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Abstract 

The two-dimensional comb lattice Ci is a natural spanning tree of the Euclidean 
lattice Z 2 . We study three related cluster growth models on C-i- internal diffusion 
limited aggregation (IDLA), in which random walkers move on the vertices of C2 until 
reaching an unoccupied site where they stop; rotor-router aggregation in which parti- 
cles perform deterministic walks, and stop when reaching a site previously unoccupied; 
and the divisible sandpile model where at each vertex there is a pile of sand, for which, 
at each step, the mass exceeding 1 is distributed equally among the neighbours. We 
describe the shape of the divisible sandpile cluster on C2, which is then used to give 
inner bounds for IDLA and rotor-router aggregation. 

Keywords: growth model, comb lattice, internal diffusion limited aggregation, rotor- 
router aggregation, divisible sandpile, asymptotic shape, random walk, rotor-router walk. 
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1 Introduction 

Let G be an infinite, locally finite and connected graph with a chosen origin o £ G. Internal 
diffusion limited aggregation (IDLA) is a random walk-based growth model, which was 
introduced by DlACONis and Fulton |DF91j . In IDLA n particles start at the origin of 
G, and each particle performs a simple random walk until it reaches a vertex which was 
not previously occupied. There the particle stops, and from now on occupies this vertex, 
and a new particle starts its journey at the origin. The resulting random set of n occupied 
sites in G is called the IDLA cluster, and will be denoted by A n . 

IDLA has received increased attention in the last years. In 1992, Lawler, Bramson and 
Griffeath [LBG92 showed that for simple random walk on Z d , with d > 2, the limiting 
shape of IDLA, when properly rescaled, is almost surely an Euclidean ball of radius 1. In 
1995, Lawler |Law95j refined this result by giving estimates on the fluctuations. Recently 
several improvements have been obtained. Asselah and Gaudilliere [AGIO! lAGllaj 
proved an upper bound of order log(ra) for the inner fluctuation 5i and of order log 2 (n) 
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for the outer fluctuation 5o in all dimensions d > 2. In [AG llbj they improve the upper 
bound on the inner fluctuation to a/ log (radius), for d > 3. Independently, and by different 
methods, Jerison, Levine and Sheffield |JLS12tlJLSllj . proved also that both Si and 
5o are of order log(n) for IDLA on Z 2 and of order \/log(radius) for d > 3 . 

Rotor-router aggregation is a deterministic version of IDLA, where particles perform rotor- 
router walks, which are deterministic analogues to random walks. They have been first 
introduced into the physics literature under the name Eulerian walks by Priezzhev, 
D.Dhar et al[PDDK96J. At each vertex of the graph G, we have an arrow (rotor) 
pointing to one of the neighbours of the vertex. A particle performing a rotor-router walk 
first changes the rotor at its current position to point to the next neighbour, in a fixed 
order chosen at the beginning, and then moves to the neighbour the rotor is now pointing 
at. In rotor-router aggregation each particle performs a rotor- router walk until it reaches 
an unoccupied site, where it stops. Then a new particle starts at the origin, without 
resetting the configuration of rotors. The resulting deterministic set R n of n occupied 
sites is called the rotor-router cluster. 

Rotor-router aggregation on the Euclidean lattice Z rf has been studied by Levine and 
Peres [LP09J, who showed that the cluster R n is a ball in the Euclidean distance. On 
the homogeneous tree Landau and Levine [LL09j proved that, provided the start con- 
figuration of rotors is acyclic, the rotor-router cluster forms a perfect ball with respect 
to the graph metric, whenever it has the right amount of particles. Kager and Levine 
[KL10J studied the shape of the rotor- router cluster on a modified two dimensional lattice, 
which they call the layered square lattice. In each of the known examples the limiting 
shape of rotor-router aggregation is the same as the one for IDLA, but with much smaller 
fluctuations compared to IDLA. 

In order to prove inner bounds for the above models, we use a third growth model, the 
so-called divisible sandpile, which has been introduced by Levine and Peres [LP09J as a 
tool for studying internal growth models on Z rf . In the divisible sandpile model each vertex 
can have an arbitrary amount of mass. If a vertex has mass at least 1, it is called unstable 
and it can topple by distributing the mass exceeding 1 equally among its neighbours. At 
each timestep a vertex is chosen and toppled if it is unstable. Provided every vertex is 
chosen infinitely often, the masses converge to a limiting distribution < 1. The set of 
vertices with limit mass equal to 1 is called the divisible sandpile cluster. If we start with 
a mass of n concentrated at the origin, the corresponding sandpile cluster will be denoted 
by S n . 

All three growth models have very similar behaviour. This was first noticed by Levine 
and Peres [LP091 ILP10] for the case when the state space is an Euclidean lattice. Com- 
puter simulations suggest that the connection between the three growth models holds in 
wide generality, but only partial results are available for other state spaces. All these 
three models have the so-called abelian property, which makes them amenable to rigorous 
analysis. In the case of IDLA and the rotor-router model this means that, if we let several 
particles run at the same time, instead of one after another, it is irrelevant for the end 
result in which order the particles make their moves. In the case of the divisible sandpile 
model, it means that the limiting distribution is independent of the order in which vertices 
topple. 

The aim of this paper is to study the three aggregation models described above on the 
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comb lattice C2, which is the spanning tree of the two-dimensional Euclidean lattice Z 2 , 
obtained by removing all horizontal edges of Z 2 except the ones on the x-axis. 

The graph C2 can also be constructed 
from a two-sided infinite path Z (the 
"backbone" of the comb), by attaching 
copies of Z (the "teeth") at every vertex 
of the backbone. 



We use the standard embedding of the 
comb into Z 2 , and use Cartesian coor- 
dinates z = (x, y) G Z 2 to denote ver- 
tices of C2. The vertex o = (0, 0) will be 
the root or the origin; see Figure [TJ For 
functions g on the vertex set of C2 we will 
often write g(x,y) instead of g(z), when 
z = (x,y). 



Figure 1: The comb C2 



While C2 is a very simple graph, it has some remarkable properties. For example, the 
so-called Einstein relation between the spectral-, walk- and fractal-dimension is violated 
on the comb, see Bertacchi |Ber06| . Peres and Krishnapur (KP04J showed that on 
C2 and other similar recurrent graphs two independent simple random walks meet only 
finitely often. Random walks on C2 have been studied by various authors, the first being 
Havlin and Weiss jWH86j and Gerl [Ger86| . 

The paper is organized as follows. In Section [2] we introduce some notations and basic 
facts which will be used through the rest of the work. Section [3] is dedicated to the study 
of the divisible sandpile on the comb Ci- We show in Theorem |3 . 5| that the sandpile cluster 
S n on C2 has up to constant fluctuations the shape 
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Section [4] deals with IDLA on €2- Using the results obtained for the sandpile model, we 
prove an inner bound for IDLA, which is of the type ([!]). Finally, in Section [5j we give an 
inner estimate for the rotor- router model on C2 which is weaker than the result obtained for 
IDLA. For a fixed initial configuration of rotors the exact shape of the rotor-router cluster 
on the comb has been obtained by the authors in [HSllJ using a purely combinatorial 
approach. 



2 Preliminaries 



Let (G, E(G)) be an infinite, undirected and connected graph, with vertex set G, equipped 
with a symmetric adjacency relation ~, which defines the set of edges E{G) (as a subset 
of G x G). We write (x,y) for the edge between the pair of neighbours x, y. In order to 
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simplify the notation, instead of writing (G, E(G)) for a graph, we shall write only G, and 
it will be clear from the context whether we are considering edges or vertices. Let o £ G 
be some fixed reference vertex called the origin. For x,y £ G, let d(x, y) be the length 
of the shortest path from x to y. Also, write d{x) for the degree of x, i.e, the number of 
neighbours of x. For a subset icGwe denote by 

dA = [x £ G \ A : By £ A with x ~ y} and diA = [x £ A : By A with x ~ y} 

the f outer) boundary respectively the inner boundary of A. 

Let P = (p(x, y)) x y( - G be the one-step transition probabilities of the simple random walk 
on G, i.e., p(x, y) = l/d(x) if y ~ x and otherwise. We write X t for the position of the 
random walker at the discrete timestep t. Probabilities will be written as P, in particular 
P x denotes the probability of a random walk which starts at x £ G. Similarly E and E x 
will denote expectations using the same convention. For y,z £ G the Green function is 
defined as 

oo 

G(y,z)=E y [Y,l { x t = z} , 

i=0 

and represents the expected number of visits to z of the random walk Xt started at y. 
For a subset A C G, write Ga for the Green function of the random walk stopped upon 
leaving the set A. That is, if r = min{i > : Xt ^ ^4}, then 

T-l 

G A (x,y)=E x [Y, t {X t =y} ■ 
For a function / : G — > K, its Laplace operator Af is defined as 

A function / : G — > M is called superharmonic on a set ^4 C G if A/ < 0, and harmonic if 
A/ = 0, for all x £ A. For a function g : G — > M, define its least superharmonic majorant 
as 

s(x) = inf {/(x) : / superharmonic , / > 5}. 
Remark that the function s is itself superharmonic on G. The following is widely known. 

Lemma 2.1 (Minimum principle.). If f is a superharmonic function on G and there exists 
x £ G such that f{x) = mine /, then f is constant. 

3 Divisible Sandpile 

Let C2 be the comb as in Figure [TJ and let /iq be an initial mass distribution on C2, i.e., a 
function hq : C2 — > M + with finite support. The divisible sandpile is a sequence (fik)k>o of 
mass distributions, which are created according to the following rule. At each time step 
k, choose a vertex x £ €2- If fJ-k(%) > 1 ; the pile of sand at x is unstable and topples, 
which means that x keeps mass 1 for itself and the remaining mass p,^ (x) — 1 is distributed 
equally among the neighbours y of x, that is, according to the transition probabilities 
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p(x,y) of the simple random walk on Given a mass distribution [i^ at time k and a 
vertex x G C2, the toppling operator can be defined as 

T x Hk(y) = Mfe(y) + »k(x)d(y)Ad' x (y), for y G C 2 , 

where <%(y) = and afc(y) = max{/ifc(y) — 1, 0}. Let (xt)k>o be a sequence of vertices 
in C2 called the toppling sequence, which contains each vertex of C2 infinitely often. Then 
the mass distribution of the sandpile after k steps is defined as 

Mfc+l = T Xk Hk = T Xk ■ ■ ■ T Xo /j,q. 

Hence, fik+i(y) is the amount of mass present at y after toppling the sites xq, ■ ■ ■ ,Xk in 
succession. One of the tools that will be used throughout this work in various incarnations 
is the so-called odometer function, which was introduced by Levine and Peres [LP09J . 

Definition 3.1. The odometer function Vk is defined as 

Vk(y)= Pjiv) - H+i(v) = Yl aj(y),yeC 2 , 

j<k:xj=y j<k:xj=y 

and represents the total mass emitted from y during the first k topplings. 

For simple random walks on C2 it is easier to work with the normalized odometer function 
Uk(x) = ^^y- Lemma 3.1 of Levine and Peres |LP09| can be easily adapted to our 
case, in order to show that, as k goes to infinity, and converge to limit functions \i 
and u respectively. Define 

S = {x G C2 : fi(x) = 1}. 

The set S is called the sandpile cluster with initial mass distribution fiQ. The limit func- 
tions fj, and u satisfy 

H(x) = fJ,o(x) + d(x)Au(x), for all x G S, (2) 

and 

fJ>(x) < I, for all x G €2- (3) 

From @ and Q it follows that 

Au(x) = — *— (l - /U (aO) for all x G 5, (4) 
d(x) 

and u{x) = 0, if x ^ S. The following result provides a method for solving this free 
boundary problem. For a proof, see once again Levine and Peres [LP091 Lemma 3.2]. 

Lemma 3.2. Consider a function 7 : C2 — > M with 

Aj(x) = (l - (Jq{x)), forallxeC 2 . (5) 

Then the normalized odometer function u of the sandpile satisfies u = 7 + s, where s is 
the least superharmonic majorant of —7. 



Lemma 3.2 gives a representation of the odometer function which is independent of the 
toppling sequence. 

Remark 3.3 (Abelian property). The limit u of the normalized odometer function and 
the sandpile cluster S are independent of the toppling sequence (xk)k>o- 
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3 DIVISIBLE SANDPILE 



3.1 Divisible Sandpile on the Comb 



With the help of Lemma 3.2, we shall next describe the limit shape of the sandpile cluster 
on the two-dimensional comb C 2 ■ Consider an initial mass distribution llq concentrated at 
the origin o, that is liq = n ■ S a , and denote by 

S n = {z G C 2 : li(z) = 1} 

the sandpile cluster, and by u n the limit of the normalized odometer function for this 
choice of initial distribution. We use another simple fact about u n ; for a proof see [LP09, 
Lemma 3.4]. 

Lemma 3.4. If x G S n \ {0} and y ~ x with d(o,y) < d(o, x), then u n (y) > u n (x) + 1. 

By Q, the normalized odometer function satisfies 

Au n (z) = -^(l - n ■ 5 (z)), for z e S n . (6) 

The odometer function u n can be reduced to odometer functions of suitable divisible 
sandpiles on Z, which are easy to compute. Let u n be the normalized odometer function 
of the divisible sandpile on Z, with initial mass distribution concentrated at 0, that is, 



l~lq = n ■ 5q. By Remark 3.3 it is clear that the sandpile cluster S n on Z in this case is a 



symmetric interval around the origin 0. In order to compute u n , by Lemma 3.2 we need 
to construct a function 7„ : Z — > WL with Laplacian given in (§. It is easy to check that 
j n defined by 

7n(y) = 2 (if I ~ \ ) ' ( 7 ) 
satisfies the required property. Since 7„ is nonnegative, the constant function is a 



superharmonic majorant of — j n . Hence, by Lemma 3.2 we have u n < j n . Now, consider 
7n : C-2 — > K with 

Jn(x,y) = Jn x (y), for (x,y) G C 2 , (8) 

where n x G K for all x G Z. The quantities n x can be interpreted as the total amount of 
mass that ends up in the copy of Z that is attached to the vertex (x, 0). Then j n satisfies 
([5]) if and only if 

n x = n- t {x=oy + % x _ x (0) - 2% x (0) + % x+1 (0) (9) 

holds for all x G Z. From ([7]) and ([9]), and using the fact that n x = n_ x by symmetry, we 
get the following recursion for the numbers n x 

n o = n +\ n \-\ n li (10) 
n x = - ^n 2 x + ^n 2 x+1 , for x > 0. (11) 



Equation (11) has strictly positive solutions as quadratic polynomials of the form 



2 n 9t 2 + 4 
n x = -x 2 -t-x-\ ^— , withiGM. (12) 
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3.1 Divisible Sandpile on the Comb 




Figure 2: Two plots of 7„ for n = 1000. The graphic on the left is superimposed 
with contour lines representing the sets B n for various values of n. In the density plot 
on the right, dark areas represent small values. By construction, the finite area which is 
surrounded by the local minima of j n coincides with the region S n covered by the sandpile. 



By the initial condition (10), the parameter t satisfies the equation 

3 



n 



which has one real root given by 



t° -\ 1, 

16 12 ' 

20. 



t = T(n)---T{n)-\ 



(13) 



with T(n) = \/2187n 2 + 125 + § n ) 3 - By a series expansion around n = oo, one 
obtains 

/ 21 4 I n 1 /^ 0(1). (14) 



Therefore, the function 7 n (x,y) = ^n x {v) with n x defined by (12) and (13) satisfies the 
conditions of Lemma \3.2\ with = n ■ 5 Q . See Figure [2] for a graphical representation of 
7„. We are now ready to prove the limit shape for the divisible sandpile on C2. 

Theorem 3.5. Let S n be the divisible sandpile cluster on C2, with = n - 5 . Then there 
exists a constant c > such that, for n > uq : 



where 



and 



Bn—c c S n cz B n -\- Cl 



(x,y)€C 2 : =ll^] 1/2 <n 1 / 3 



(15) 



3 \ 2/3 

2 



1 / 3 



2 V 2 



1/3 
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Proof. The upper bound S n C B n + C : The mass distributions n x are nonnegative for all 
x, therefore 7 n is nonnegative, and this implies that the constant function is a super- 
harmonic majorant of — j n . Thus, by Lemma |3.2[ j n is an upper bound of the odometer 



function u n . Moreover, Lemma 3.4 implies that u n decreases by a fixed amount on the 
sandpile cluster S n when we move away from the origin. Therefore, in order to get an 
upper bound for S n , it suffices to calculate the minima of 7 n along each infinite ray starting 
at o = (0, 0). By the symmetry of C2, it is sufficient to consider only the first quadrant. 

Consider the rays which lies entirely on the positive x-axis. We have 7 n (x, 0) = The 
minimum of this function is attained at x m i n = 

jt, with t given in (|13|). Using the series 



expansion ( 14 ) of t we get 

2/3 



X 1 ! 



kn 1/3 + <D(l), with *;=(-) , (16) 



which is also an upper bound of S n on the x-axis by Lemma 3.4, since 7n(|_27ninJ > 0) is 



bounded by a constant which is independent of n, and smaller than 1/10. 

To calculate the extent of the sandpile cluster on the "teeth", we need to compute the 
minima of 7 ra in the y-direction. On each "tooth" of the comb, j n is a quadratic polynomial 



which attains its minimum at y m i n (x) = Moreover, 7(x, [jAmnMJ) < 1/2. Using (12) 
and a series expansion around infinity we get 

y mi n(x) = I (V/3 - |) 2 + \x - InV 8 - ^xn^ 3 + 0(1), 

1 /3 

where I = 2 (|) • By the estimate in the x-direction we know that (x,y) 6 S n only if x < 
^min- Thus, using the expansion (16) for x m i n we obtain (x,y) G S n if \x\ < kn 1 ^ 3 + 0(1) 



and \y\ < I (n 1 / 3 - f) +0(1), for n > no- This proves the upper bound S'n C i3 n _)_ c . 

T/ie lower bound B n - C C S 1 ,!.' On each infinite ray the minimum of y n (z) is smaller than 
a constant a > 0, independent of n. Also, from the upper bound, we have u n (z) = for 
all z G dB n+c . Hence it n (,z) — 7n(^) > — o for ai l 2 S dB n+c . Since the function n n — 7„ is 
superharmonic, by the Minimum Principle, it attains its minimum on the boundary and 
the inequality u n (z) — y n (z) > — 0, holds for all z G S n+C . Thus 7„ — a is also a lower 
bound of the odometer function on B n+C , which gives the inner estimate B n - C C S n , for 
some constant c. □ 

The next corollary follows directly from the proof of the theorem. 

Corollary 3.6. Let u n be the normalized odometer function of the divisible sandpile on 



C2, with initial mass distribution /io = n ■ 5 , and B n C C2 defined as in (15). Then there 
exists a constant < a < 2, such that, for all n> n$ and all z G C2 

(j n (z) - a)l Bn < u n (z) < j n (z). 



4 Internal Diffusion Limited Aggregation 

Let (JQ) . gN be a sequence of independent and identically distributed simple random walks 
on the comb C2, with common starting point Xq = o. Then X\ represents the position of 
the i-th particle at time t. 
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Definition 4.1. Internal diffusion limited aggregation (IDLA) is a stochastic process of 
increasing subsets (Aj) igN of C%, which are defined recursively as A\ = {o} and for i>2 

F[Ai = Au{x}\ Ai_ x = A] =P[X % ai =x], 
where a 1 = inf{t > : X\ ^ Aj_i} is the first exit time of the random walk X\ from Aj_i. 



The IDLA cluster is build up one site at a time. 
That is, suppose that we already have the cluster 
Aj_i after % — 1 particles stopped, and we want 
to get Aj. For this, the i-th particle X\ starts at 
o, and evolves as long it stays inside the IDLA- 
cluster Aj_i. When X\ leaves Aj_i for the first 
time, it stops, and the point outside of the cluster 
that is visited by X\ is added to the new cluster 
Aj. The set Aj is called the IDLA-cluster of i 
particles. 

Figure [3] shows IDLA clusters on C 2 with 100, 500 
and 1000 particles. 

We will prove the following shape theorem for 
IDLA on C 2 . 



Figure 3: IDLA cluster 

Theorem 4.2. Let A n be the IDLA cluster of n particles on C 2 . Then, for all e > 0, we 
have with probability 1 



£> n (i_ e ) C A n , for all sufficiently large n, 



where 



(17) 



1/2 



<n 1 /3 



and 



The set B n is the same as the limit shape of the divisible sandpile from Theorem 3.5 The 



proof of Theorem 4.2 



Levine and Peres 



uses ideas of Lawler, Bramson and Griffeath |LBG92j and of 
LP10J. Following [LBG92J, we introduce the stopping times 



r* = min{t >Q:X\$ B n } and t\ = min{f >§:X\ = z\, 
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4 IDLA 



for z G B n . Consider the probability that a fixed vertex z G B n does not belong to the 
IDLA cluster A n , which can be written in terms of the stopping times defined above as 



a < t z 

i<n 



Hence, by the Borel-Cantelli Lemma, convergence of the series 



£ Yl n^A n \, (is) 

n>n zeB n(1 _ e) 



is a sufficient condition for Theorem |4.2[ Fix now n and z G B n and consider the random 
variables 



N = J2t {T i <ai} , M = J2t {rl<K} and L = ^1 {(T , 



<T i <T i \ 



Then A represents the number of particles that visit z before leaving the cluster. The 
variable M counts the number of particles that visit z before leaving B n , and L is the 
number of particles that visit z after leaving the cluster Ai but before leaving B n . Remark 
that if L < M, then z G C2 belongs to A n . Moreover N > M — L. Therefore, in order to 
estimate F[z A n ], we just need an upper bound for F[M = L\. For any fixed number a 

F[z £ A n ] = F[N = 0] < F[M - L = 0] < F[M < a or L > a] 

< F[M < a] +F[L > a}. ^ 

We shall show that for a suitable choice of a, the probabilities F[M < a] and F[L > a] are 



small enough, such that the series (18) converges. The derivation of a suitable value of 
a will be done differently, not like in the case of Euclidean lattices, studied by Lawler, 
Bramson and Griffeath in [LBG92J, who used classical asymptotics for the Green 
function stopped on a ball. Since in our case, the Green function Gg n stopped on B n is not 
directly available, we use the odometer function of the divisible sandpile as a replacement, 
as suggested by Levine and Peres in [LP10] . For simplicity of notation, we shall write 
G n (y, z) instead of Gg n (y, z). The random variable M is a sum of i.i.d. indicator variables, 
with 

E[M] = nF [r z < r n ) = n^4- (20) 

^nKZ, Z) 

Even though L is a sum of dependent indicator variables, following [LBG92J, L can be 
bounded by a sum of independent indicators as follows. Only those particles with X 1 1 G B n 
contribute to L and for each y G B n there is at most one index i with X l i = y. The 
corresponding post-r^ random walks are independent. In order to avoid dependencies in 
L, enlarge the index set to all of B n and define 

l y 

{r 2 <r„}' 



yeB n 

where the indicators l v correspond to independent random walks starting at y. Then 
L < L, and the expectation of L is given by 

E[L] = £ F y [r z < r n ] = £ G n (y,z). (21) 

yeB n ' y&B n 
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Now (19) can be rewritten as 

F[z i An] < F[M <a}+ F[L > a}. (22) 

We shall relate the random variables L and M with the odometer function of the divisible 
sandpile. For this, consider the function f n : dB n U B n — > R, 

f n (z) = 9lM m -L]. (23) 
Set h x (y) = Gn }u^ , then for x,y £ B n one has 



d(y) 

1 



if x = y 

d( x y \ (24) 

0, ifx^y 
By linearity of the Laplace operator, f n solves the following Dirichlet problem 

Af n (z) = -^-(l-n-5 {z)), fOTZEBn 
a{z) 

f n (z) = 0, for z £ dB n - 

Recall that the odometer function u n of the divisible sandpile with initial distribution 
Ho = n - 5 solves the same Dirichlet problem ^ on the sandpile cluster S n (whose shape 
is given by B n by Theorem 3.5). 

The uniqueness of the solution of a Dirichlet problem gives that f n = u n on the set B n , 
and u n is approximated (up to an additive constant) by the function 7„ defined in Q. 
Since u n > 0, it follows that f n (z) > 0, for all z G B n , which is equivalent to E[M] > E[L]. 

We will use the following large deviations estimate for sums of independent indicators. 
For a proof, see Alon and Spencer |AS921 Cor. A. 1.14]. 

Lemma 4.3. If N is a sum of finitely many independent indicator random variables, then 
for all A > 0, 

F[\N - EN\ > XEN] < 2e' c ^ EN , 
where c\ is a constant depending only on A. 



In order to find an upper bound for the right hand-side of ( 22 ) we use the previous Lemma 
and choose A > and a such that 

(1 + A)E[Z] < a < (1 - A)E[M]. (25) 

Hence A has to satisfy the relation 

, E[M — L] Uz) , , 

< A < J- = ^-i, 26 

E[M + L] g n {z) 

and g n defined as 

g n (z)= C ^^E[M + L\. (27) 

To obtain A, we have to bound f n {z) / g n (z) away from 0. For this, we first have to calculate 
the function g n which is, like f n , defined on dB n U B n . 
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4 IDLA 



4.1 The function g n 



By (20) and (21), the function g n is the solution of the Dirichlet Problem 

n ■ 6 (zj), 



&9n{z) 



1 



d(z. 



("1 



9n{z) = 0, 



for z £ B n , 
for z £ dB n , 



(28) 



and can therefore be obtained by solving some linear recursions. For simplicity, we first 
shift the set B n by kn 1 / 3 in the direction of the positive x-axis. This shifted set will be 
denoted by B^, which is the set of all (x, y) G C2 with < x < 2kn 1 ^ 3 and 



< 



for < x < kn 1/3 , 



\y\ < 



(2kn 



1/3 



for kn 1 / 3 < x < 2/cn 1 / 3 . 



On the shifted set we define the function g n : B n — > R, by 

9n( x , V) = 9n{x + kn 1/3 , y), 



(29) 



which solves the same Dirichlet problem ( |28[ ) on B l n with the origin moved to (kn l l 3 ,ti). 
By symmetry of g n , it is enough to compute for vertices (x,y) with < x < kn 1 / 3 and 
y > 0. For z = (x,y) £ B l n , with y 7^ 0, the Laplace Ag!^(z) is equal to —1/2, hence on 
each "tooth" of the comb, g l n satisfies the linear recursion 

25n(^> V) = 9l,(x, y + 1) + gl (x, y - 1) + 1, 
which has the general solution 



9n(x,y) = 77 (y - y 2 ) + c 1 (x) + yc 2 (x), 



(30) 



where c\{x) and c 2 (x) are functions of x, to be determined. For (x, 0), (x, 1) G C 2 , we have 

SnO^ ) = ci(x) and #*(x, 1) = ci(x) + c 2 (x). (31) 

From (28) we have the boundary conditions g l n (0,0) = and g i n (2kn 1 ^ 3 , 0) = and for 
< x < kn 1 / 3 , we have g^(x,x 2 /3) = 0. On the other hand, from equation (30), we get 



^(x,x 2 /3) 



6 



x 



+ ci^) + y c 2(^) = 0, 



which implies that the function c 2 (x) can be written as 



c 2 (x) 



1/x 2 



1 



,ci [X . 



2 V 3 

Moreover, on the x-axis the Laplace operator of g l n satisfies 

-i, if x / kn 1 ! 3 

-j(n + l), if x = fcn 1 / 3 . 



A^(x,0) 



(32) 



(33) 
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4.2 IDLA inner bound 



For x 7^ fcn 1 / 3 , that is, when (x, 0) is not the center of B^, we have 

<£(x + 1,0) = 4^(z, 0) - tfjx - 1, 0) - 2<^(x, 1) - 1, 



and using (31) we obtain 

c\(x + 1) = 2c\{x) — c\{x — 1) — 2c 2 (x) — 1, 



which together with (32) gives an equation for ci, namely 



ci(x + 1) = (2 + ^)ci(x) - d{x - 1) 



X 

x* 1 J ' ' 3 
This has an explicit solution as a polynomial of degree 4, given by 



c\{x) 



1 

18 



x 4 + fcc 3 



1 

36' 



-x 



(34) 



where 6 is a free parameter which can be computed using the other boundary conditions 
for g*. Since Ag^(/cra 1//3 , 0) = — |(n+ 1), using equations (31), (32), and (34), we obtain 



5K + 27n 
18(1 + 3K 2 )' 

3 A 2/3 



4.2 



Since we 



where K = kn 1 ^, and the constant k = (| ) is the same as in Theorem 
are interested in the form of for n sufficiently large, we expand b around n = 00, giving 

b{n) = -n l / z + 0(n- 1 / 3 ). 

Putting everything together we get g n (x,y) = g^^kn 1 ^ 3 — \x\, \y\), with 

<7n(z, y) = (^« 1/3 + 0(n" 1/3 )) (x 3 - 3xy) + ^(3y - 18y 2 - 2x 4 - x 2 + 12x 2 y). 



4.2 IDLA inner bound 

We are now able to conclude the proof of Theorem |4.2| 

Lemma 4.4. For all e > there exists n £ , such that for all n > n £ and all z £ £> n (i„ e ) 

e < E[M - L] 



4 E[M + L] 



Proof. By (26), one needs to study the function X n (x,y) 



We have 



(\v\-%y 



2ci(x) + (2c 2 (x) + l)y-y 2 ' 



where ci(x), c 2 (x) and n x are defined in (p4|), p2[) and (12), respectively. It suffices to 



consider the first quadrant. For every fixed x, the function A n (x,y) is decreasing in y for 
< y < From the proof of Theorem 3.5 we already know that 



2 



(35) 
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4 IDLA 



For < e < 1 consider the set 



(x,y) G C 2 : |x 



< (1 - e)^ 1 / 3 and \y\ < (1 - e)Z (n 1 / 3 - ^ | 



Obviously £> ni£ C B n for all e, hence & is well defined on this set. Furthermore, by (35), 



Li 



is also decreasing on B n ^ £ as a function of y, for all e > and n big enough. This means 
that it is enough to study ^ at the inner boundary of £> n ,e- For this, let z = (x, y) G B Ut£ 

with \y\ = (1 — e)l (n 1 / 3 — J?H be such a boundary point. Then 



lim 



/n(*) 



e e 
> 



n ^ 9n{z) 4-e 4 
The statement follows from the fact that for each e > one can find an e' > such that 

^n(l-e) C B n , £ '. 



□ 



Proof of Theorem 4-2 Recall that we need to show the convergence of the series (18). Fix 
z G S n n_ e ). We set A = | > in Lemma 



4.4 



and choose 



o = (l + A)E[L] = (l + i jE[L] 



in equation (25). Apply now Lemma 4.3 to M and L. Recall also that E[M] > E[L]. Then 



[M < a] + F[L > a] = P M < ( 1 + - )E[L] 



< 4exp{-c A E[L]} < 4exp|-c A 



L> (1 + JE[L] 

9n(z) ~ f n {z) 



G n (z,z) 

where c\ is a constant depending only on A. Hence, for all n > n e , we have 

9n(z) ~ f n {z) 



E E p[^^„]<4yj yj e Xp {- 

n>n E 2eB n (!_ e ) n>n E zeB n ( 1 _ e ) 



-ca- 



G n (z,z) 



}■ 



(36) 



In order to estimate the stopped Green function G n (z, z) upon exiting B n , with z = (x, y), 
note that 

\y\<b n (x) : =/(nV3_M) 2 . 

We have the trivial upper bound G n (z, z) < 2GA(y, y) where Ga is the Green function of 
the simple random walk on the integer line, stopped at the interval A = [ — b n (x), b n (x)j . 
Using Proposition 1.6.3 and Theorem 1.6.4 from Lawler |Law91j . this can be bounded 
by 

n , v b n (x) 2 -y 2 ( 1/3 |a;|\2 

For every e > 0, the function g n (z) — f n (z) is again decreasing on every non-crossing path 
which starts at o and stays inside B n n £ y Hence, it attains its minimum on the inner 
boundary diB n n £ \ of B n (i_ £ ). Taking limits, we get for every sequence z n = (x,y n ) with 
x fixed and z n G diB n n £ \ 

,. 9n{Zn) — fn{Z n ) k, . 

hm w = - 2 - e e, 
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4.3 The recurrent potential kernel 



and for the sequence z' n = (x n , 0) with x n = kn l /^(l — e) 1 / 3 

Um 9n{z' n ) - Uz'n) = k. _ ^ _ _ _ 1/3 . 

n— >oo n I 4 

Hence for all e > and n big enough 

min ( 5n (z) - / n 0)) > C £ • n 4/3 , 



for a constant C £ which depends only on e. Since, by (37) the stopped Green function 
Ga(z, z) is of order 0(n 2 / 3 ), this implies 

9n(z) — fn{z) „, 2/3 

mm — — — — > C £ • n 1 . 

2eB n( i_ e) & n (z,z) 



Hence, (36) can be bounded by 

Yl F i z( £ A n\<4j2 nexp{-c x C' £ n 2 ^} < oo, 

n>n e zEB n n— e ) i%>n e 
which concludes the proof. □ 



4.3 The recurrent potential kernel 

In the recurrent lattice case Z 2 , the limiting shape of IDLA is derived using estimates for 
the recurrent potential kernel which is defined as follows 

n 

A(x, o) = lim V (F [X t = o] - P x [X t = o}), for x G C 2 . 

t=0 

See Lawler, Bramson and Griffeath [LBG92j for more details. For some constant 
N > 0, the level sets of the potential kernel are sets of the form {x G C2 : A(x, o) > N}, 
and the level sets of the Green function are of the form {x G C2 : G(x, o) > iV}. 

In all previously known cases, the limiting shape of the IDLA-cluster is determined by 
the level sets of the potential kernel in the recurrent case and by level sets of the Green 
function in the transient case. Nevertheless, this is not the case for comb lattices C2, even 
if the simple random walk is recurrent. In order to show this, let us consider the generating 
function of the potential kernel 

00 

A(x, o\z) = {Wo[X t = 0}- P x [X t = o\y 

t=0 

= G(o,o\z) — G(x,o\z), 

where G(x,y\z) = X^o^I^* = y\ zt ls ^ ne generating function of the Green function of 
simple random walk on Ci- Using standard techniques for generating functions one gets 
for x = (xi, X2) G C2 

G(x,o\z) = F 1 (z)^F 2 (z)^G(o,o\z) 
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5 ROTOR-ROUTER AGGREGATION 



with 




1- VT^2 



y/2y/l - z 2 + y/l 



— z- 



2 



Z 



z 



and 



G(o, o\z) 



V2 



a/1 - z 2 + VI 



— z 



2 



See Bertacchi and Zucca |BZ03j for details. Therefore 



A(x,o\z) = G{o,o\z) (l - F 1 (z)^F 2 (z)\ x ^ 



and lim 2 _ s>1 - A(x,o\z) = 2|xi|. By Abel's theorem for power series, the potential kernel 
on the comb C2, if it exists, has to be equal to A(x, o) = 2|xi|. Hence we have an example 
of an IDLA cluster whose behaviour is not given by the level sets of the potential kernel. 



5 Rotor-Router Aggregation 

A rotor-router walk on a graph G is defined as follows. For each vertex a; fix a cyclic 
ordering c(x) of its neighbours, i.e., c(x) = (xo, x±, . . . , xm x )-i), where x ~ x% for all 
i = 0, 1, . . . , d(x) — 1. The ordering c(x) is called the rotor sequence of x. A rotor 
configuration is a function p : G — > G, with p{x) ~ x, for all x £ G. Hence p assigns to 
every vertex one of its neighbours. A particle configuration is a function a : G — > No, with 
finite support. If a(x) = m > 0, we say that there are m particles at vertex x. A particle 
located at a vertex x with current rotor p(x) = Xj, performs a rotor-router walk like this: 
it first sets p(x) = X{+i, where addition is modulo d(x), and then it moves to Xj+i. 

Rotor-router aggregation is a deterministic process of increasing subsets (Riji^ of G de- 
fined recursively as R\ = {o}, and 



where Z{ is the first vertex outside of Ri—\ that is visited by a particle performing a rotor- 
router walk, started at o. The particle stops at Zi, and a new particle starts its tour at 
the origin, but without resetting the rotor configuration. The set R n of the occupied sites 
in G is called the rotor-router cluster of n particles. The odometer function ur{x) of 
rotor-router aggregation is defined as the total number of particles which are sent out by 
the vertex x during the creation of the rotor-router cluster R n . 

In this section we study rotor- router aggregation on C2 , and we give an inner bound for the 
cluster R n which holds for arbitrary initial configuration of rotors and is independent of the 
rotor sequence. The approach below relies on an idea of Holroyd and Propp [HP10], 
who use rotor weights in order to prove a variety of inequalities concerning rotor-walks 
and random walks. 



Ri = Ri-i U {z^ for i > 2, 
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5.1 Rotor Weights 



5.1 Rotor Weights 

Let G be a locally finite and connected graph, gq : G — > N an initial particle configuration 
with finite support, and po : G — > G an initial rotor configuration with po(x) = xq for 
all x 6 G, that is, all initial rotors point to the first element in the rotor sequence c(x). 
Routing particles in the system, such that at each time step t exactly one particle makes 
one step of a rotor- router walk, gives rise to a sequence (pt, at)teN °f r otor- and particle- 
configurations. To each of the possible states (pt, at) of the system, we will assign a weight. 

Fix a function h : G — > R. Define the particle weights at time t to be 

W P (i) = 2<7((i)Mi), (38) 

and the rotor weights of single vertices x € G as 

, ,. (O, for/c = 

= < (39) 
- 1) + h(x) - h(x kmodd(x) ), for k > 0, 

where x» is the i-th neighbour of x in the rotor sequence c(x). Notice that, for k > d(x), 

w(x, k) =w(x,k — d(x)) - d(x)Ah(x). (40) 

The total rotor weights at time t are given by 

W R (i) = ^ u>(x,u t (x)), 

where is the odometer function of this process, that is, the number of particles sent 
out by x in the first t steps. Note that po is chosen in such a way that, if i = ut(x) mod d(x), 
then Xj = pt{x) for all t > and x G G. 

It is easy to check that the sum of particle- and rotor- weights are invariant under routing 
of particles, i.e., for all times t,t' > 

W P (t) + W R (t) = W P (0 + W R (0- (41) 

For rotor-router aggregation on G we start with n particles at the origin, that is, do = n-5 , 
and we route a particle only if there is at least one other particle at the same position. The 
process terminates when no two particles are at the same position. Denote by t* = t*(n) 
the number of steps it takes to finish the process, and by (o~t*,pt*) the final configuration. 
By the abelian property, the configuration (at* , pt* ) does not depend on the order the 
particles made their steps, and by definition at*(x) = l| xe ^ n }. 

We use the following weight function: for some y 6 G, let h y : G — > R given by 

h(x) = h y (x) = (42) 

where G n is the Green function of the simple random walk on G, stopped upon exiting 
the sandpile cluster S n with initial mass distribution po = n • 5 Q on G. Take now some 
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5 ROTOR-ROUTER AGGREGATION 



y E S n . Here S n is the sandpile cluster of the divisible sandpile on some general graph 
G. Recall that the Laplace of h y (x) on G is given by (24). The particle weights at the 
beginning are 

W F (0) = nh y (o), (43) 

while the rotor weights are Wr(0) = 0. At the end of the process, i.e., at time t* when 
the rotor-router cluster R n is formed, we have 



w P (i*) = m*) ^ E 



(44) 



since h y is equal to outside of S n . For the rotor weights we get from (40) 



w R (0 = E 



d(a?) 



d(x)Ah y (x)) + w(x,k x ), 



(45) 



where u/j is the rotor odometer function and fc^ = ur(x) mod d(a:). By (24) and (39) 



w R (0 



ur{v) 

d(y) 



+ EE(^ 

x£R n i=0 



hy(Xi)) 



(46) 



+EEiw-w 



Hence by the invariance of the total weights ( |41[ ), we obtain 

£ (ndo(x) - l)^(x) < ^ + £ £|M*0 " ( 47 ) 

Denote by u(y) the lefthand side of ( |47[ ). Then v(y) solves the Dirichlet problem 

(Av(y) = -^(l-nS (y)), for y E S n , 
\v(y) = 0, for y S n . 

By Q , the normalized odometer function u n of the divisible sandpile on G with initial mass 
distribution fiQ = n ■ 5 Q satisfies exactly the same Dirichlet problem, hence v(y) = u n (y) 
and we get the following result, which compares the odometer function ur of rotor- router 
aggregation with the odometer function u n of the divisible sandpile. This result holds for 
any locally finite and connected graph G. 

Proposition 5.1. Let u n be the normalized odometer function of the divisible sandpile 
with initial mass distribution fj,Q = n • S 0) and ur the odometer function of rotor-router 
aggregation with n particles starting at the origin o £ G. Then, for all y EG, 



n„(y)<^#+W R (y) 



with 



Wr(!/) = E E 



d(y) 

G n (y,x) G n (y,z) 



x£S„ z~x 



d{x) 



d(z) 



(48) 
(49) 
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5.2 Rotor-router Aggregation on the Comb 



Levine and Peres |LP09] derived an inequality similar to (48) in the case of Z using 



a different method. For trees, Wr(j/) can be expressed in terms of the expected distance 
from the starting point of a random walk to the point where it first exits S n . 

Proposition 5.2. If G is a tree and d(-,-) is the graph distance on G, then 

W n (x) = 2E x [d(x,X T )] -2, 
where T = inf ji > : Xt S n }, and (Xt) is the simple random walk on G. 



Proof. For y ~ z let N yz be the number of transitions from y to z before the random walk 
exits S n . Then 

E]N _ N , G n (x,y) _ G n (x,z) 
x[ yz zyJ " d(y) d{z) • 

See also |LP} Proposition 2.2] for more details. Since G is a tree, the net number of 
crossings of each edge is smaller or equal to one, i.e., 

\& X [N yz -N zy \\< 1. 

We consider G as a tree rooted at x, and denote by tt x>z the shortest path from x to z. 
For y / i, write y~ for the parent of y, i.e., the unique neighbour of y that lies on the 
shortest path ir x ,y With this notation we get 



E 

y,z£S„ 
y~z 



G n (x,y) G n (x,z) 



d{y) 



d(z) 



2 Ex[JV, 



-JV - 
y y iy yy v 



y,zeS n 
y~z 



S/G5„ 
y^x 



where the last equality is due to the antisymmetry of N yz — N zy . Let 

C y = {z £ S n : y G ir XtZ } 

be the cone of y. The random variable N y - y — N yy - is either zero or one, the latter if the 
random walk exits S n in the cone C y , hence 



z . 



y&Sn 

y^x 



y&Sn zeCy 

y^x 



For all z G dS n we have #{y G S n \ {x} : z G C^} = d(x, z) — 1, therefore 

W R (x) = 2 F X [X T = z] (d(x, z)-l)= 2E X [d(x, X T )} - 2, 

zddS n 

which completes the proof. 



□ 



5.2 Rotor-router Aggregation on the Comb 



Since C2 is a tree, by Proposition 5.1 and 5.2, one needs an upper bound for the expected 
distance from the starting point of a random walk (Xt) on C2 to the point where it first 
exits S n , in order to derive an inner estimate of the rotor-router cluster. Recall that on 
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5 ROTOR-ROUTER AGGREGATION 



C 2 , the sandpile cluster S n has the shape B n given in Theorem 3^ Using the trivial upper 
estimate 

K z [d(z, Xt)} < max {d(z, w) :w£ dS n } = \x\ + |y| + ln 2/3 , (50) 

-1 in 

with z = (x, y) and I = h (| ) as in Theorem 3.5 we can show the following inner bound. 



Theorem 5.3. Let R n be the rotor-router cluster of n particles on C 2 . Then, for n > uq 
and for any initial rotor configuration and choice of rotor sequence, we have 



B n C Rni 



where 



B, 



• n = { (x, y)eC 2 : \x\ < kn 1 ' 3 - c x n x l\ \y\ < I (n 1 / 3 - |)' 



+ c 2 x - c 3 n 



where 



2/3 



I) ' 



2 V 2 



1/3 



and c\, C2 and C3 are constants. 



Proof. By the definition of rotor-router aggregation, jz 6 C 2 '■ ur(z) > 0} C R n , and by 
Proposition 5.1 together with Proposition 5.2, we have for vertices z = (x,y) 



d(z) 



> u n {z) - 2E Z [d(z, X T )] + 2> u n {z) - 2(|x| + |y| + In 2 / 3 ) + 2, 



where the last inequality is due to (50). By Corollary 3.6, we have a lower bound of the 
sandpile odometer u n for z S S n 

ln(z) - a < u n (z), 

where a is a positive constant smaller than 2, and j n is the function defined in ([8]). Thus, 
to derive an inner bound, it suffices to check for which z = (x, y) £ S n the inequality 



7„(x,y)-2(|a;| + |y|+/n 2 / 3 ) >0 



(51) 



holds. By symmetry it is enough to consider x,y > 0. We first check inequality (51) on a 
"tooth" of the comb, that is, for a fixed x. The function j n is given as 



7n(x,y) 



1 / n a 
V 



2 V 2 

where n x is the amount of mass that ends up in the x-"tooth" of the sandpile. Since x 



is fixed, we can treat constant. Hence the right hand side of (\51h is a quadratic 

polynomial in y with smallest root 



Vx 



2 + 



n, 



4 + yre 2 / 3 + 2n x + Ax. 



Substituting n x as calculated in (12), an expansion around n = 00 gives 

- 2 2 + V& 



y.i 



In 2 ' 3 



1 

27 



n 1/3 x + — + 



1 /3 X 

-x + c\n > ' — c 2 



n 
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5.2 Rotor-router Aggregation on the Comb 



Since (x,y) G S n , we have the bound x < kn 1 / 3 , hence 

I (n 1 / 3 



y.v 



x\ 2 , 2 + VE 1/3 



(52) 



for n > no, and a positive constant c. To get a bound on the x-axis, we calculate for which 
x > the inequality y x > is satisfied. Since y x is a polynomial of degree 2 in x this is 
easy to do, and again by series expansion around n = oo we obtain 



< k ■ n 1 / 3 - cgn 1 / 6 , 



for n > uq. The inner bound for R n now follows from (53) together with (52). 



(53) 



□ 



Figure [4] shows the inner estimate of the rotor-router 
cluster from Theorem 5.3 in comparison to sandpile 
cluster S n , for n = 1000. The white area is the area where 



the inequality (51) is valid, and corresponds to the set B n 
of Theorem 5.3 The colouring is based on the value of 



the right-hand side of (51) 



The inner bound could be improved if one has a sub- 
stantially better upper bound for E 2 [d(z, Xt)\- For reg- 
ular graphs, one can also give an universal inner estimate 
for rotor-router aggregation, which relates the rotor-router 
cluster to a divisible sandpile cluster with a smaller mass. 
Using the methods of Levine and Peres in |LP09j one 
can deduce the following. 

Proposition 5.4. Let G be a regular graph with degree d 
and root o, and let R n be the rotor-router cluster of n par- 
ticles starting at o. Further, let S n be the divisible sandpile 
cluster with mass distribution ^o(x) = n ■ 5 (x). Then 
S n /(2d-i) c R n . 




-10-5 5 10 



Figure 4: Inner bound 



Remarks. In the IDLA model, for an upper bound of the type A n C B n (i +£ ^, we still do 
not have sufficient tools. In all previous studied cases, the IDLA cluster grows uniformly, 
and this makes easy the study of random walks. This is of course violated in our case, 
since the set B n in ([!]) grows with rate n 1 / 3 in the x-direction, and with rate n 2//3 in the 
y-direction. The harmonic measure for random walks stopped when exiting the set J3 n , 
was studied in [HSllj . by making use of a special rotor-router type process. This is the 
first example where IDLA aggregate is not a set of uniform harmonic measure. In [HSllj, 
we also give subsets of C2 with uniform harmonic measure. 



Recently Duminil-Copin et al. |DLYYllj developed a method for proving an outer 
bound for IDLA without needing a harmonic measure estimate, provided an inner bound 
is known. Nevertheless, their method cannot be applied in our setting, since one of the 
required regularity conditions (weaker lower bound) is not satisfied on the comb. 
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